Moisture induced Ageing in Granular Media 
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We present experiments showing that a granular system of small beads does 
exhibit ageing properties : its maximum stability angle is measured to increase 
logarithmically with resting time, i.e., the time elapsed before performing the 
measure. We show that humidity is the crucial ingredient responsible for this 
behaviour : while ageing effects are important at intermediate humidity, they 
disappear at vanishing humidity. On the basis of these experimental results, we 
propose a model based on the activated condensation of liquid bridges between 
the beads. Within this picture, we are able to reproduce both the waiting-time 
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q \ and humidity dependence of the ageing properties. 
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^ ■ In his pioneering treatise published in 1773 [T|, Charles de Coulomb pointed out the 

peculiar properties of granular systems. His decisive remark was to recognize the links 
existing between the statics of granular systems and the friction properties of two solids at 
contact : the equilibrium and stability of granular edifices could be discussed in terms of 
the friction between the different part or layers of the system ||. In particular, he deduced 
a relationship between the rest angle of the pile, 8q, and the static friction coefficient fi s 
of the medium, in the form tan#o — Us- One knows by now that this relation is only a 
good first approximation and that a wide variety of behaviours can be observed in granular 
systems ||. The same is true for dry friction problems : while the subject remained poorly 
understood during more than 150 years, a lot of new results emerged in the last four decades. 



In particular, one of the intriguing behaviour of friction is the existence of ageing properties 
of the static friction coefficient, fi s : the latter is found to increase logarithmically with 
resting time, i.e., with the time during which the system was left at rest before pulling it (in 
order e.g. to probe fi s ). Now, if Coulomb's analogy between solid friction and the statics of 
granular systems is further pursued, a question arises : can ageing be observed in granular 
systems too ? 

In order to verify this assertion, we have studied the effect of waiting time on the angle 
of first avalanche 9 W of a granular system of small (typically 200/im) spherical glass beads 
contained in a rotating drum, following a protocole described in fig. l.a. Astonishingly, 
and in close analogy with dry friction, a logarithmic ageing of the maximum static angle 
9 W is observed, as shown in fig. 2. a. In other words, granular media exhibit ageing. The 
dependence on bead size has been studied and ageing was not observed for beads with a 
diameter larger than 0.5 mm, except for very large humidities. On the other hand it is known 
that humidity is an important parameter in the description of granular media [Q]. Addition 
of small quantities of wetting liquid has been shown to change enormously the repose angle 
of a pile ; a brief discussion of moisture effects can even be found in Coulomb's treatise 
|J. We thus repeated the experiments for various humidities P v /P sat . As shown on figure 
2.b, humidity appears to be the crucial parameter which controls the ageing of 9 W : no 
ageing is observed for small humidity, while the "strength" a of the ageing effect increases 
in an enormous way with humidity. 

This humidity dependence leads to the intuitive idea that the ageing effects originate from 
the condensation of small liquid bridges between the beads. Indeed, liquid bridges induce 
a significant cohesion force between the beads, which can increase the friction between 
different layers of the granular system and result in a higher value of 8 W . However, the 
physical justification of this intuitive behaviour is not obvious at all. Let us first consider 
the idealized situation of two smooth beads in contact : the capillary adhesion force exerted 
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by a small liquid bridge connecting the beads is given by the product of the liquid surface 
tension times the radius of the beads [] : ~ 2ii'jR, whatever the size of the (small) bridge 
||. Therefore, the adhesion force has no humidity dependence, and would be moreover able 
to stick all beads together 0]. This apparent difficulty disappears if one takes into account 
the roughness of the beads. At the liquid-vapour equilibrium, the radius of curvature r eg of 
the liquid meniscus is fixed by the so-called Kelvin relation || 

— = pi k B T log ^* = pi Ap (1) 
where pi is the density of the liquid and the ratio defines humidity. Under ambient 

*sat 

conditions, this yields a nanometric order of magnitude for r eq ! A crucial consequence is 
that liquid bridges are able to form only in nanometric interstices. Because the beads 
are not smooth at the nanometer scale, the wetted region does not spread over the whole 
possible area it would occupy if the beads were smooth, and the cohesive force is reduced in 
proportion. Another counter-intuitive behaviour is the slow evolution in time of the cohesive 
properties, measured through the time dependency of 9 W . This indicates that condensation 
of liquid bridges takes place over very long time scales. This feature is in agreement with 
experiments using Surface Force Apparatus (SFA) and AFM techniques 0, which show that 
an interstice between two solid surfaces can remain for a very long time in a metastable 
dry state while the equilibrium state would be a condensed liquid bridge. 

On the basis of all these considerations, we propose a model based on the condensation 
of interstitial liquid bridges via an activated process. In this letter, we shall only outline 
the main features of the physical mechanism and a more detailed presentation will be pub- 
lished elsewhere by the authors. Let us consider two "rough" beads at contact. Under the 
assumption of an activated process, the time needed in order to condense a liquid bridge in 



x The following expression assumes that solids do not deform because of the capillary force. If 
solid deformation occurs, the numerical prefactor changes slightly 0]. 
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an interstitial volume is r ~ r exp (j^f)) with r a microscopic time and AE an activation 
energy barrier. Nucleation occurs preferentially at the level of a nano-asperity, and the acti- 
vation energy is accordingly AE ~ A/i pia^e, where Ap = p sat — p g — A;_bT log (P sa t/ 'P v ), pi 
is the density of the liquid, e the gap between the surfaces at the level of the nucleating site 
and <2q a typical nucleation area (see fig. l.b). Now, since both beads are rough, one expects 
the nucleating sites to exhibit a broad statistics of gaps e between solid surfaces and the 
activation times are accordingly widely distributed. After a given time t w , only the bridges 
with an activation time r ac t smaller than t w have condensed. These were therefore formed at 
the nucleating sites with a gap e verifying e < e max (t w ) = ^f-^^r log (^)- Once a liquid 
bridge has condensed, it fills locally the volume surrounding the nucleating site, until the 
Kelvin equilibrium condition for the radius of curvature, eq. (0), is met. Thus, because of 
roughness, only a part of the total wettable area is indeed wetted at a given time t w . The 
corresponding fraction f(t w ), is proportionnal to the number of activated bridges, yielding 
in first approximation f(t w ) ~ e max (t w ) / X, with A the typical width of the distribution of 
distances between the surfaces. The capillary adhesion force is thus reduced by the same 
factor as compared to the perfectly smooth case, leading to 

* 7 d r-Lj log(^) (2) 

log ^ \T J 

where d = 2tcR/ '(Xpia^) is a distance taking into account the geometrical characteristics of 
the contact. Now, by reproducing Coulomb's argument for the stability of the surface layer 
in the presence of this additional adhesive force, one obtains the following implicit equation 
for 9 w (t w ) : 

a(P v ) ft w \ 

tan9 w (t w ) ~ tan0 o H j- log — J (3) 

cos 8 w (t w ) \t J 

with a(P v ) = cto/log-^f 1 and ao is a numerical constant depending on the characteristics 
of the beads. Thus, by plotting tan9 w (t w ) as a function of log (t w )/ cos9 w (t w ), one should 
obtain a straight line : as shown on fig. 2. a, this expectation is indeed in very good agreement 



with experimental results. Moreover, as exhibited on fig. 2.b., the increase of the slope, 
a(P v ), of this line with humidity is in good agreement with the theoretical prediction eq. 
(0). Our model is thus able to reproduce both the waiting-time and humidity dependence 
of the measured ageing properties of a granular system. 

Our results highlight the crucial role of humidity in the statics of granular systems. 
We have shown that the latter do exhibit ageing properties in an analogous way to those 
encountered in dry friction for the static friction coefficient /i s . The analogy is in fact even 
more striking, since a similar effect of humidity has been reported in friction between rocks 
||, though not extensively studied : in ref. ||, the "standard" ageing properties of fi s 
were found to disappear at vanishing humidity, just like in the case of a granular system ! 
Similar behavior were even observed in indentation experiments too |J. It would be thus 
interesting to check to which extent the analogy between dry friction and granular media 
is indeed pertinent. The hope is to propose an alternative unsder standing for the ageing 
properties observed in solid friction. 
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Figure l.a : Description of the experimental setup. The glass beads fill 25 % of the volume 
of a cylinder (diameter 100 mm, thickness 13 mm) which can rotate around its axes. Two 
different systems have been extensively studied, a "polydisperse" one, 140fim < d < 260/zm 
and a "monodisperse" one 200/im < d < 250fim, d being the diameter of the beads. The 
walls of the cylinder are made of glass, and the shape of the bead's pile is recorded with a 
video camera. Experiments are performed at room temperature, under controlled humidity 
(defined as the ratio between the vapour and saturated pressure of water, P v /P sat ). Before 
starting experiments, the system is prepared by rotating the drum during approximatively 
12 hours. Then the ageing properties are investigated for various values of P v /P sa t using 
the following protocol : (i) first, the system is put in motion for a few turns (typically 
three); (ii) the end of this "motion period" defines the origin of waiting time t w = 0, after 
which the system is left at rest; (iii) after a given time (ranging from 10 to 10 4 seconds), a 
slow rotational motion (with velocity less than one turn per minute) is transmitted to the 
cylindrical drum. The angle 9 W at which the first avalanche takes place is measured from the 
slope of the beads just before the avalanche. The same procedure is repeated for different 
waiting times t w and a whole curve 9 W vs. waiting time t w can be constructed. 

Figure l.b : Schematic drawing of the contact at the nanometer scale between two 
micrometric asperities on the beads. Inside most of the contact region, the solid surfaces do 
not really touch eachother at the molecular scale, and capillary condensation occurs in the 
interstice left between the surfaces. 

Figure 2. a : Logarithmic ageing of the angle of first avalanche. As explained in the text 
after eq. fl3|), the ageing property is best analysed by plotting tan 9 w (t function of 

log 10 (t TO )/ cos9 w (t w ), for different values of the water vapor pressure P v . From bottom to 
top, humidity is : 15% (triangles), 27% (pentagons), 36.1% (squares) and 45.5% (circles). 
In these measurements, times runs typically over the range t = 5 s up to t = 5000 s. The 
dotted lines are least-square fits of the experimental data, whose slope is identified with 
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a(P v ). 



Figure 2.b : Variation of the slope a(P v ) characterizing the ageing behavior of the 
first avalanche angle (see text eq. (|3|)) with humidity P v /P sat . Open dots correspond to 
the polydisperse system, and filled dots to the monodisperse one. The dashed line is the 
theoretical prediction a = a / log(P* at /P v ), where a = 0.079 and P* at = 0.68 P sat . The 
validity of the previous theoretical prediction is best confirmed by the measured linear 
dependence of 1/a as a function of \og(P sat / P v ) (see inset) : in this plot, a linear least- 
square fit thus provides unambiguously both values of ao and P* at . The lowering observed 
in the saturating pressure P* at might be an effect of the long-ranged attraction forces exerted 
by the walls, and/or of dissoluted species in condensed water. 
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Figure la. Moisture induced Ageing in Granular Media 
Corresponding author : L. Bocquet 
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Figure lb. Moisture induced Ageing in Granular Media 
Corresponding author : L. Bocquet 
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Figure 2a. Moisture induced Ageing in Granular Media 
Corresponding author : L. Bocquet 
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Figure 2b. Moisture induced Ageing in Granular Media 
Corresponding author : L. Bocquet 
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